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DIGITAL COMPUTER 



LIBRARY ROUTE© t 7 ~ 18B 



TITLE 

TYPE 
DURATION 



DESCRIPTION 



METHOD OF USE 



Iterative Estimation of CcnBunalities % Principal Axes 
Method of Factor Analysis 
Complete program 

Dependent on number of iterations and order of matrix. 
Each iteration takes one minute for order 11 and approx- 
imately five minutes for order 18 (the largest order 
matrix for which this program can "be used). 
The code makes a principal axes factoring of the given 
matrix, revest imates the commonalities, substitutes the 
new estimates for the old, and prepares for another 
principal axes factoring. The latent roots and new com- 
monality estimates are punched during each iteration » 

th 
The i eommunality estimates are 



(i) h 2 
J 



(1)J2 



"J* 



J = 1, 



s P 



■where p is the order of the matrix, r if the hypothesized 
rank, and the factors ^"'f '. of the i iteration have "been 
ordered, in terms of their contribution to variance. 

1. The program tape is input until it stops. 

2. The problem tape is input on the black switch, 

3. The program tape is placed in the reader just after 
the word last read by the computer. Start with black 
switch to Stop Disable * Iteration will begin* When 
it is desired to end run, raise the black switch to 
Obey * The computer will stop on a "2^ GF9 W order 

at the end of the iteration it is going through. 
Now start with the white switch* The lower off-* 
diagonal and diagonal elements of the correlation 
matrix will be punched by rows with the last eom- 
munality estimates as diagonal entries* This latter 
output will simplify preparation of the problem tape 
that will be required for further iteration. Ee*o~ 
gram ends on an n 0F" order. 



PREPARATION OF TAPES 



COMMENTS 



EXAMPLE 



-2- 



The program tape can be copied from the master tape . 
The problem tape has the parameters followed by two 
copies of the lower off -diagonal and diagonal elements 
of the matrix punched by rows . The character "N" is 
punched after the last element of each copy. The 
elements of the matrix are punched in the form required 
for Library Routine N 3« 

The parameters at the beginning of the problem tape 
are punched as follows? 

d J p F r L 
where d is the number of digits desired In the printed 
latent roots and communality estimates t j> is the order 
of the matrix, and r is the rank of the matrix for which 
communality estimates are wanted. 

The squared multiple correlations may be suggested as 
initial wsmmunality estimates. 

The elements of the matrix should be. multiplied 
(scaled) by one-tenth or ~by one-hundredth so that the sum 
of squares over the entire matrix is less than one-half. 

The latent roots and communality estimates will be 
scaled by the same constant as the elements of the matrix 
(either one-tenth or one -hundredth). 
We would like communality estimates for a set of four 
■variables for which we have the following correlation 
matrix i 



1.00 


•73 


„H 


.25 


•73 


1.00 


.k 9 


.36 


M 


•*9 


1.00 


-31 


.25 


.36 


•37 


1.00 



We decide to retain two factors. We scale the elements 
of the matrix by one-tenth since the sum of squares for 



^ 



he formula is 



R 11 -1 

Hii—Wli i 



th 



R 



li 



where R J is the element in the i row and j 



.th 



column of the inverse of the correlation matrix^ the inverse can be 
obtained with Library Routine M 2 . 



-3- 



the matrix will then he less than one-half, Because we 
want cesHBHtmalities to four decimal places we Must print 
five places since all answers will he scaled hy the con* 
stant_, one~tenth. 

The prohlem tape is punched as follows? 



5 

J 
k 
F 
2 
L 

+0543 

+ 073 
+ Qhh 
+ 025 
spaces 
+05^3 
+ 073 
+ okh 

+ 025 



results printed to five places 



natrix is of order four 



assumed rank is two 



+0585 
+ 0+9 
+ 036 



+0585 
+ 0-^-9 

+ 036 



+0298 
+ 037 



+oi8iir 



+0298 
+ 037 



+0181H 

The results for the first two iterations will he punched as? 



-01952 
+05882 



+18103 
+06775 



-01352 
+03555 



+01271 
+01892 



-91^59 +18721 +OI615 -0077^ 
+05993 +07231 +03650 +OI847 



DATE £ 

CODED BY 



l*j 195? 



SSg 



APPROVED 




LOCATION 


ORDER 


JTOT& PAGE 1 ; 

- — - ■- . . ! 


Routine XI - 


18 


Decimal Order Input 




00 165K 











00 F 

00 285F 

00 171K 




First storage location 







80 F 
00 F 




(-1) 


- Constants 


1 


00 IF 










L5 ( )F 


by 50F 


Test -word 






00 100K 








I 


loutine H 3 - 


23 


Decimal Number Input 




00 20K 




Input routine 





*H 2F 




i 




iH F 


from 8L 


1 


1 


81 kF 
LO 28L 


from 5L 


i 
1 
\ 


2 


32 5L 




i 




iA 28L 








3 


.50 F 
7^ 28L 








k 


S5 F 
40 F 






— Input independent parameters 


5 


26 U. 

^2 6l 


from 26 






6 


L5 F 

40 (00 )F 


by 5L 






7 


F5 2F 
kO 2F 








8 


LO 2% 
32 OL 








9 
• 


50 kF 
75 ^F 








f 10 . 

i 


S5 F 

L4 4F 








i u 

-i 


10 IF 
Lk 165F 









K7 



LOOTIOH 



13 
14 

15 
16 

17 
18 

19 
20 
21 
22 

23 
24 

25 

26 

27 
28 



ORDER 



HOTES 



PAGE 2 



K 7 



4o 7F 
26 30L 
L5 4f 
40 162F 
40 15F 
LO 5F 
40 166F 
L5 5F 
40 i64f 
L5 3E 
40 161F 
00 20F 
46 218F 
46 226F 
50 4f 
75 4f 
S5 F 
L4 163F 
4o 167F 
36 42F 
00 20F 
46 25L 
L5 4F 
00 IF 
4o 5F 
92 567F 

40 ( )F 

50 25L 
26 100F 
26 80F 
00 F 
00 F 
00 F 
00 10F 



from 31L 



Waste 



from 92F 
by 22L 



— Form and store dependent parameters 



Punch delay characters 



— Input correlation matrix 



LOCATION 


ORDER 


IT0TE3 Eft.SE 3 


29 


80 F 
00 3F 






- Constants for input routine 


30 


k2 1J2F 
ho 163F 


from 12L 






31 


26 13L 
00 F 










00 258K 






Routine P 2 - 


32 from 227F 
from 219F 
from 276F 


Print (A) "with or -without Sign to n 
Places 




00 80K 




Checls: input of correlation matrix 





L5 167F 
42 2L 


from 46f 






1 


LO I65F 
L4 163F 








2 


22 12L 
L5 (00 )F 


from 13L, 6L 






3 


l4 59F 
^4-0 59F 


"by OL, 4L 




-Form sum. of matrix just input 


4 


F5 2L 
42 2L 








5 


LO l4L 

32 6l 








6 


22 2L 











L5 59F' 


from 5L 


i 


7 


LO 60F 




} 


8 


40 F 
L3 F 




p Compare sums from the tm> correlation 
! matrices 




34 999F 




_i Return to input if suras do not agree 


9 


L5 59F 




~i 




40 60F 








10 


41 59F 
92 451F 






Print stroke character 


11 
s 


L5 15L 
40 11L 






" Prepare to input correlation matrix 
once again 


12 


26 45F 




— ' 






42 10F 


from 2L 


h 



K 7 



LOCATION 


OKDER 


HCTES PAGE 4 


13 


42 i4l 




—Addition to comparison sub-code 




22 2L 






14 


22 12L 








L5 (00 )F 


"by 13L 


— Constants 


15 


OF F 








26 45F 




- 




00 233K 




Transfer correlations > 





L5 163F 
42 % 


from 276F 




1 


L5 165F 








42 4L 




-Plant addresses 


2 


L5 167F 
00 20F 






3 


46 kL 








36 kL 


Waste 




4 


L5 ( ;f 


from 8L,"by JL " 






kO ( )F 


lay 1L, 6L 




5 


19 18F 

F4 4L 






6 


kO kL 
k6 9L 




—Transfer correlation 


7 


LO 9L 
36 243F 






8 


26 4L 








36 9L 


^feste 




9 


L5 ( )F 


by 6L 


Test word 




kO ( )F 


by OL 






00 24jK 









L5 163F 
4a % 


from 230F 




1 


42 13L 
4l I69F 






2 


4i i68f 


from 12L 






L5 168F 


from 9L 




3 


LO 169F 






i 


40 170F 







K 7 



LOCATION 


ORDER 


NOTES PAGE 5 


4 


L3 170F 
32 12L 


1 




5 


32 5L 


¥a£te 






41 ( )F 


"by OL, 6L 




6 


F5 5L 
42 5L 


from 14l 




7 


42 13L 
F5 i68f 




— Generate identity matrix 


8 


4o l68F 
LO l62F 






9 


36 10L 
22 2L 






10 


F5 169F 
4o 169F 


from 9L 




11 


LO 162F 
36 213* 






12 


26 2L 








L5 171F 


from 4l 




13 


50 IF 








40 ( )F 


"by 1L, 6l 




14 


26 6l 








00 F 








00 222K 




Put together and print latent roots 





L5 165F 
42 3L 






1 


00 20F 
46 3L 




— Prepare addresses 


2 


41 168F 








92 135F 




Panch tw> line*f eed characters 


3 


L5 ( )F 
40 ( )F 


from 10L 
"by 6L, 9L 

by OEL, 10L 


Transfer a latent root 


4 


50 (8)F 


"by 38F 






50 4L 




— Print the latent root 


5 


26 258F 








F5 168F 


" 


] 



K7 



Locmair 


OKDER 


HOTES PAGE 6 


6 


40 168F 
LO l62F 






7 


36 175F 
F5 168F 






8 


00 20F 
L4 3L 




-Increase comcters and test for end 


9 


46 3L 
F5 3L 






10 


42 3L 








26 3L 








00 175K 




Set (p-k) smallest roots eqral to zero 





4l 168F 


from 229F 






41 173F 




-dear counters 


1 


41 I69F 








L5 165F 


from 13L 




2 


L4 169F 






3 


42 6l 
42 7L 
42 9L # 




-Plant address of next root 


4 


42 16L 
50 1023F 






5 


L5 I65F 








L4 173F 




-Plant address of smallest root so far 


6 


42 8l 


._ 






L3 ( )F 


hy 2L 


1 


T 


32 11L 








LI ( )F 


fey %x 




8 


32 15L 






9 


L4 ( )F 

36 14L 


by 6L 


-Compare enrrent root 




14 ( )F 


&y 31 




10 


4o OF 
L3 OF 






11 


36 14L 




_ 




F5 169F 


from 71* 151 
17L 





K7 



LOCATION 



ORDER 



NOTES 



PAGE 7 



KT 



12 

13 

Ik 

15 
16 

17 
18 

19 

20 
21 



from %* 11L 



from 8L 



by ^L 



from 13L 



1+2 169F 
LO l62F 
32 17L 
22 1L 
L5 169F 
1+2 173F 

L5 168F 

1+2 168F 

iH ( )F 

22 11L 

L5 165F 

L*+ 173F I 

k2 19L j 

92 IF j 

kl ( )F ! by 18L 

F5 168F 

k2 168F 

LO 166F 

32 22L 

00 197K 

22 175F 

in 168F 
92 13IF 
19 38F 
1+0 170F 
L5 167F 
LO 170F 
k2 20L 
L5 20L 
F^ 168F 
k2 20L 
1+1 170F 
1+1 I69F 
L5 165F 



from I96F 



from 24l 



-Increase and test coutnter 



Go back to finish a round of comparing 



j- Current root is smallest root, so far 
! this round 



-One of the (p-k) smallest roots has 
been f 012nd 



Estimate, transfer, and print commonalities 



Punch one linefeed character 
2~ 59 to 170F 



-Pre-set an address 



LOCATION 


ORDER 


1 — — : 1 

HOKES PAGE 8 K 


7 


42 12L 


j \— Set addresses 




50 162F 


1 




8 


75 168F 
L5 163F 


i 




9 


S4 F 
42 10L 








10 


42 11L 
50 (00 )F 


by 9L, 16L 
from 19L 






11 


32 11L 


Waste 








7J (00 )F 


hy 10, 16L 




— Form a squared factor loading and add to 


12 


40 1023F 
50 (00 )F 


by JL, 15L 




the sum across a row of the factor matrix 


13 


7J 1023F 
L4 170F 








l4 


40 170F 

F5 12L 








15 


42 12L 
F5 10L 








16 


42 10L 








17 


42 11L 
F5 169F 
4o 169F 






- Prepare to finish a commonality estimate 


18 


LO 162F 
32 19L 




- 




19 


22 10L 










00 IF 


from 18L 




20 


L5 170F 
40 (00 )F 


by 3L, 5L 






21 


50 (8)F 
50 21L 


by 38F 




— Print and store the commonality estimate 


22 


26 258F 
F5 168F 








23 


40 168F 
LO 162F 






— Test for end 





LOCAIIOH 


OKESl 




ffOTES PAGE 9 


2k 


36 276F 








26 in. 


Return to estimate next coaHonality 




00 800K 






Routine X 2 - 


108 Shifting Sua Cheek 




F68HL6S695 


Siam check mnriber 




2k 20N 






Problem tape 


Ls inserted here and read| the Program tape is then 




re-ins^rted at this point. 






00 276K 









92 to 


from 22XF Punch delay characters 




? 4- 255F 


Stop before starting new iteration 


1 


L5 167F 
M2 2L 


h 


I 


2 


92 191F 








L5 (00 )F 


from 6L, by 1L^ 




3 


50 7F 

50 3L 


j 




k 


26 258F 




1 




F5 2L 


Print correlation matrix with last 


5 


te 2L 
LO 8L 




communality estimates 


6 


36 7L 

22 2L 






7 


92 770F 
OF F 


from 6L 




8 


92 191F 








L5 SK 








00 20K 






] 


Routine M k *« 


136 1 


Closed Eigenvalues and Eigenvectors 




00 152K 






3 


Roirtin: S 1 - 


116 


Sqxaare Root Rout in- 




26 ^33H 


. 





K7 



